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networks 

E.C. Behrman, R.E.F. Bonde, J.E. Steck, and J.E Behrman 



Abstract — Entanglement of a quantum system depends upon 
relative phase in complicated ways, which no single measure- 
ment can reflect. Because of this, entanglement witnesses are 
necessarily limited in applicability and/or utility. We propose 
here a solution to the problem using quantum neural networks. 
A quantum system contains the information of its entanglement; 
thus, if we are clever, we can extract that information efficiently. 
As proof of concept, we show how this can be done for the 
case of pure states of a two-qubit system, using an entanglement 
indicator corrected for the anomalous phase oscillation. Both the 
entanglement indicator and the phase correction are calculated 
by the quantum system itself acting as a neural network. 



I. Introduction 

The scientific process asks questions of the universe. Much 
of the education of a scientist consists of learning, painstak- 
ingly, how to ask those questions. If the question is ill-posed, 
the answer will be ambiguous or even unintelligible; even 
if it is well-posed, the question needs to be the one whose 
answer is pertinent to whatever problem it is we wish to 
solve. Simulations are one good route to finding out about 
the clarity and worth of our questions, since variables can be 
easily controlled and extraneousness determined. 

Quantum imformation is the systematic use of the fun- 
damental quantum mechanical nature of the universe to do 
calculations that are either very difficult or even impossible to 
do with a classical computer. Exactly what this allows us to 
do, and how we can make use of that nature, turns out to be 
a very difficult problem. The quantum system itself contains 
that information. But how do we exploit those abilities? How 
do we even find out what they are? 

In this paper we use a Quantum Neural Network (QNN), 
which necessarily complex-valued, to calculate the essential 
attributes of the quantum information. While this is done in 
simulation, the method, as outlined, could be implemented 
experimentally to create a quantum computer, programmed 
through neural network training, to calculate its own entan- 
glement and phase. 
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II. Entanglement 

Determination of entanglement is one very good example. 
This is an important question because entanglement is what 
allows us to do quantum computations (TJ; it is an essentially 
quantum question since no classical system can have quantum 
correlations; thus, it is an excellent testbed for this general 
problem. The quantum system which is our quantum computer 
"knows" what its entanglement is, but how do we ask the 
question? And how do we know if we have the right answer? 

Given the density matrix, entanglement can actually be cal- 
culated for the two-qubit system. For definiteness, we choose 
the "entanglement of formation" Ep [2], which is the number 
of pure singlets necessary to create a given entangled state. Of 
course, we may not even know the density matrix. We could 
determine the density matrix through quantum tomography 0, 
but the number of measurements necessary goes like 2 2Ar , 
where N is the number of qubits. Thus, as the system grows 
this method becomes prohibitive. Is there instead a relatively 
cheap way to determine the entanglement? 

A number of researchers (see, e.g., (U, 0) have devised 
single measurement "entanglement witnesses", to indicate 
the presence or absence of entanglement. An entanglement 
witness, W, is a Hermitian operator such that tr(Wp) > 
for a density matrix p representing a fully separable (i.e., 
unentangled) state. Most of them are measures of "closeness" 
to a single, entangled state, and fail for other kinds of 
entangled states. This is easiest to understand by considering 
the so-called "Bell" states: 
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|*-((9 4 )> = ^[|01)-e^|10) 

For {Oi} = 0, this is a commonly used basis for the two-qubit 
system, the "Bell basis" (it is an orthonormal set). Each of 
these states is fully entangled, Ep = 1 (i.e., each could be 
created using exactly one singlet), independent of the angle 

A witness can be devised that measures closeness to any 
one of these states. For example, the Mintert witness (4) is 
given by W M = -4.tr (p <g> pV), where V = p[ 1] <g> (P.i 2) - 
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Fig. 1. The Mintert Witness, W M , set for closeness to |^-(<9 4 = 0)>, 
tested on each of the Bell basis states as functions of phase offset 0. Both 
|$+) and |<E>_) are calculated to have zero entanglement, while |^+) and 
|^_) oscillate, showing entanglement for half the range of 6. Note that the 
correct answer is that each state is fully entangled for all 6 values (i.e., the 
witness should be negative everywhere.) This oscillation is what we are calling 
"anomalous oscillation" in the entanglement witness. 



(2) 

i+ ), P- is the projector onto the antisymmetric states, and 
P + is the projector onto the symmetric states. As proposed 
this witness is set up to work well for states close to the 
singlet (antisymmetric) state, |\£_(04 = 0)), so it does not 
detect entanglement in either |<£ + ) or |$_). Worse, the witness 
increasingly fails to detect entanglement in |\£_(#4)) as #4 
increases. See Figure [T] which shows Wm ( set for closeness 
to the EPR singlet, |^_(0 4 = 0))) , calculated for each of the 
Bell states as functions of their respective 6. Entanglement is 
indicated by negativity of the result. The straight line showing 
the constant result of zero shows that this witness predicts 
no entanglement for either |3>+) or |3>_). The two oscillatory 
curves correspond to the and \S&-) states. At 6 = tt/2 

each becomes the other; thus, at that point, is predicted 
to lose entanglement and to become entangled. 

Witnesses can be "reset" so that they will work for other 
members of the Bell basis; for this witness, all that is necessary 
is to change the projection operators. See Figure [2| which 
shows results of this procedure: all the predictions now lie atop 
each other (as they must do by symmetry.) So, if you know 
which entangled state your unknown state is close to, you can 
use a witness which will correctly predict its entanglement. 
However this creates an obvious dilemma: you must have 
knowledge of the state before you do the measurement to gain 
knowledge of the state. 

All single measurement witnesses must show this anoma- 
lous oscillation [6], shown in Figures [T] and [2] Clearly, it would 
be of value to devise a corrector to this oscillation, so that 




Fig. 2. The Mintert Witness, Wm, set for closeness to each of the Bell 
basis states (at 9 = 0) in turn, and tested as functions of phase offset 6. All 
show exactly the same oscillatory behavior. Again, the correct answer is that 
each state is fully entangled for all 6 values (the witness should be negative 
everywhere.) 

the witness would then correctly predict that all the states of 
Equation 1 are fully entangled. With such a method, we could 
thereby determine whether a general, unknown, input state was 
entangled, and perhaps even the degree of entanglement. We 
propose here such a method, using a quantum neural network. 

III. Quantum neural network 

In previous work we have shown that a quantum system can 
be trained to deal with both parts of the problem: that it can 
be trained to output an entanglement indicator O, l8lL l9l or 
a phase indicator (7). Briefly, the details are as follows. 

For the 2-qubit system, we write for the Hamiltonian: 

H = K A CF xA + K B CF x b + £A°zA + ^B<^zB + C°zA°zB (2) 

where {a} are the Pauli operators corresponding to each of the 
two qubits, A and B, Ka and Kb are the tunneling amplitudes, 
Sa and Sb are the biases, and ( the qubit-qubit coupling. 
(Generalization to an N-qubit system is possible 151, f9l but 
will not be considered here.) The time evolution of the system 
is given by the Schrodinger equation: 

where p is the density matrix and H is the Hamiltonian. 
The parameters {K, e, (} determine the time evolution of the 
system in the sense that, if one or more of them is changed, 
the way a given state will evolve in time will also change. 
This is the basis for using our quantum system as a neural 
network. The role of the "weights" of the network is played 
by the parameters of the Hamiltonian, {K, e, (}, all of which 



JOURNAL OF JATeX CLASS FILES, VOL. 6, NO. 1, JANUARY 2007 



3 



we take to be experimentally adjustable as functions of time 
(see, e.g., |10], for the case of SQuID charge qubits.) By 
adjusting the parameters using a neural learning algorithm we 
can train the system to evolve in time to a set of chosen target 
outputs at the final time tf, in response to a corresponding 
(one-to-one) set of given inputs. Because the time evolution is 
quantum mechanical (and, we assume, coherent), a quantum 
mechanical function can be mapped to an observable of the 
system's final state, a measurement made at the final time 
t f . The time evolution of the quantum system is calculated by 
integrating the Schrodinger equation numerically in MATLAB 
Simulink, using ODE4 (Runge-Kutta), with a fixed integration 
step size of 0.05 ns ITT1 . The system was initialized in each 
input state in the training set, in turn, then allowed to evolve for 
190 ns. A measurement is then made at the final time; this is 
the "output" of the network. An error, target — output, is cal- 
culated, and the parameters are adjusted slightly to reduce the 
error. This is repeated for each (input, target) pair multiple 
times until the calculation converges on parameters that work 
well for the entire training set. Complete details, including a 
derivation of the quantum dynamic learning paradigm using 
backpropagation lfT2l in time fl3lL are given in (6) and in Q. 

A. QNN entanglement indicator 

For the entanglement indicator, we chose as our output 
the qubit-qubit correlation function fill evaluated at the final 
time, (o- Z A(tf)cr Z B(tf)) 2 , and trained the indicator using a 
training set of four (one fully entangled state, two product 
(unentangled) states, and one partially entangled state.) That 
is, the input,output) pairs are: 

input =|^(0)) (4) 
output = (o- z A(tf)cr zB (tf)) 2 — » target 

with prepared input states at zero time, and corresponding 
targets, given in Table [I] which also shows the trained val- 
ues and the entanglement of formation, calculated using the 
analytic formula for comparison. Note that the QNN indi- 
cator systematically underestimates Ep for partially entangled 
states; this is because we found through simulation that the 
net naturally trained to the target value of 0.44. See (6) for 
details. That is, we seek here not exact agreement with E F 
(in which case we would train the state -^(|00) + 1 01) + |10)) 
to a target value of 0.55), but a robust and internally self 
consistent measure, which we would hope would track well 
with an analytic measure like Ep- 

Note that the results in Table U are retrained from our 
previous work [6], this time using continuously varying func- 
tions rather than functions piecewise constant in time; training 
was, as might be expected, much more rapid, and the trained 
parameter functions more symmetric: Ka and Kb lie right on 
top of each other, as do ca and 6b- Both e and ( were simple 
oscillatory functions, of only a single frequency; however, K 
seems to exhibit two frequencies. We curvefit each function to 
a Fourier series. The agreement was very good. Coefficients 
for the fits are shown in Table [n] The trained functions, with 
their curvefits, are shown in Figures [3] and [4] 



Entanglement Parameter Curvefits 
0.15i f== 




50 100 150 200 

Time (ns) 



Fig. 3. The functions 6 a = £b and £, as functions of time, as trained for 
the entanglement indicator, and plotted with a single frequency Fourier fit. 
Each was started out (pre-training values) as a constant function: 6 a = £b — 
C = 10- 4 GHz. 



Entanglement Parameter Curvefit 
1 .956 1 1 1 . — 




50 100 150 200 

Time (ns) 



Fig. 4. The function Ka = Kb as a function of time, as trained for 
the entanglement indicator, and plotted with a two-frequency Fourier fit. K 
was started out (pre-training values) as a constant function: Ka = Kb = 
1.875 x 10- 3 GHz. 
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TABLE I 

Training data for QNN entanglement witness. 



Input state ^(0)> 


Target 


Trained 


E F 


^(|oo> + |ii» 


1.0 


0.998 


1.0 


|(|00> + |01> + |10> + |11» 


0.0 


1.2 x 10- 5 


0.0 


_^(0.5|10> + |11» 


0.0 


1.8 x 10- 4 


0.0 


^(|00> + |01> + |10» 


0.44 


0.44 


0.55 


RMS 

Epochs 




4.4 x 10~ e 
100 





TABLE II 

CURVEFIT COEFFICIENTS FOR PARAMETER FUNCTIONS K, €, f , FOR QNN 
ENTANGLEMENT WITNESS. 



fit) = a + 


aicos(ujt) + bisin 


(bJt) + CL2COsi2ojt) 


+ b2sin(2ojt) 


coefficient 


K 


epsilon 


zeta 


ao 


0.0019495 


1.014 x 10~ 4 


1.012 x 10~ 4 


ai 


-1.002 x 10- 6 


2.824 x 10- 5 


1.109 x 10- 5 


bi 


6.868 x 10- 6 


9.577 x 10- 6 


-3.96 x 10- 5 


a 2 


2.981 x 10~ 6 






b 2 


-4.562 x 10- 7 






UJ 


0.01645 


0.02674 


0.05282 


RMS 


1.069 x 10- y 


1.88 x 10~ b 


7.982 x 10-° 



This indicator gives good results for large classes of input 
states, including both pure and mixed states, and is not 
restricted to states "close" to any particular state. We have 
also 0, extended our work to the 3-, 4-, and 5-qubit 
cases, and found that as the size of the system grows, the 
amount of additional training necessary diminishes; thus, our 
method may be very practical for use on large computational 
systems. But as it is a single measurement, we too get 
anomalous oscillation. See Figure [5] which shows that the 
QNN indicator tracks very well with Ep for states of the 
form aoo|00) + ane^|ll) as a function of aoo, but predicts, 
incorrectly, that the entanglement is a function of <fi. Another 
way of looking at this is shown in Figure [6] and Figure [7] If 
we confine our testing to states with only real coefficients, the 
QNN gives an excellent approximation to the entanglement 
of formation. Figure [6] shows a comparison of the QNN 
entanglement indicator for 50,000 randomly generated states 
with real coefficients for the 2-qubit system. Agreement is 
excellent (45° yellow line is ideal); however, if we relax the 
restriction, and include complex coefficients (or, equivalently, 
nonzero phase offsets), our agreement becomes quite bad (see 
Figure [7]) So, we need to know the phase offset(s) of the 
system before we can calculate the entanglement. 

B. Phase indicator 

Since the phase offset is yet another quality of the system, 
the system should, itself, be able to supply us with information 
about its own phase. And so it can. We choose the "charge" 
basis, in which a general pure state can be written 

|*(0)> = a 00 |00) + aoie*|01) + a 10 e^|10) (5) 

+ane^|ll) 




Fig. 5. QNN entanglement indicator for states of the form aoo 1 00) + 
an e*^ 1 11), where aoo and an are both positive real, as a function of both 
relative magnitude and <f>, and compared with the entanglement of formation. 



2-qubit Pure State QNN - Real Coefficients 




Fig. 6. QNN entanglement for 50,000 randomly generated pure states of 
the form aoo|00) + aoi 1 01) + aio|10) + an |11), where aoo, a>oi, aio,and 
an are all real, as a function of the entanglement of formation. Points lying 
along the dashed yellow line are states for which the entanglement predicted 
by the QNN witness exactly matches the entanglement of formation. 



where normalization requires that 

\J a oo + a oi + a lo + a li = 1 (6) 
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Fig. 7. As in Figure [6] but with complex coefficients. 

TABLE III 

CURVEFIT COEFFICIENTS FOR PARAMETER FUNCTIONS K, 6, £, FOR QNN 
PHASE INDICATOR. 



f(t) = ao + a\cos{ujt) + b\sin(ujt) 



coefficient 


K 


epsilon 


zeta 


ao 


0.002512 


8.945 x 1CT 5 


7.445 x 1CT 4 


ai 


5.156 x 10- 5 


-1.005 x 10- 5 


-6.346 x 1CT 4 


bi 


-3.781 x 10- 6 


8.4e - 005 


1.359 x 10- 4 


UJ 


0.0658 


0.03454 


0.06402 


RMS 


7.556 x 10~ e 


3.869 x lCT b 


4.468 x lCT b 



Since an overall phase is physically meaningless we may take 
out any overall phase factor; that is, without loss of generality 
we may take the coefficient of the 1 00) basis state to be real. 
We then write each of the other coefficients as its magnitude 
times a phase factor; thus, each a nm will be a real number, 
and the phase factor, if any, will be written in explicitly. There 
are, thus, three independent phases for the two-qubit system. 
Our training set consisted of (input, output) pairs of the form 

input = |tf(0)) = 4=(l°°> + (7) 
v2 

output = \(U\V(tf))\ 2 H> target = cos 2 (>/2) 

with 1 1 different values of <p from — it to it. Training was rapid 
and good. Figure [8] shows the trained parameter functions for 
the phase indicator. Ka and Kb lie right on top of each other, 
as do €a and e b • Each was started out (pre-training values) as 
a constant function: Ka = Kb = 2.5 x 10 ~ 3 GHz, and ca = 
€b = C — 10 _4 GHz. Each trained rapidly to what seemed 
obviously to be a single oscillatory function. We therefore 
curvefit each function to a single term in the Fourier series. 
Coefficients for the fits are shown in Table [TTU 



Phase Parameter Curvefits 

3i 1 . 




50 100 150 200 

Time (ns) 



Fig. 8. The functions Ka = Kb, £a — e B, an d C as functions of time, as 
trained for the phase offsets, and plotted with curvefits to single oscillatory 
functions, as explained in the text. On the scale of the graph the curvefits are 
indistinguishable from the plotted points of the trained values. 

Agreement is good even with the fitted functions in place of 
the trained ones; and it is certainly easier to use a fixed analytic 
function than a large collection of discrete data samples. For 
testing on more general states than those of Equation 7, 
though, we needed to expand the target function. 

Because quantum computational power comes from en- 
tanglement, it is not surprising that the ability of the net 
to determine a relative phase between two basis states de- 
pends on there being entanglement between those states, and, 
indeed, the size of the signal diminishes with decreasing 
entanglement. For the Bell states with unequal magnitude, 
&oo|00) + an 1 11) , we found that if we used the target 
function 2(^ — «oo) 2a ii +2aoo&ii cos 2 (0/2) we got excellent 
results. Note that 2aoo&n, the coefficient in front of the 
cos 2 (0/2) term, is the concurrence for this state - a monotonic 
measure, like E F , of the entanglement [2 ]. Testing on states of 
the form aoo|00)+aoi|01)+ane^|ll) was equally good when 
we again adjusted for the diminished amount of entanglement 
in the target function. 

In fact, the parameter functions {Ka = Kb,ca — £bX} 
we found through training only on the eleven states of 
Equation 7 gave us all three of the phases {0, with 
symmetrically adjusted outputs and target functions (e.g., to 
find we used the output |(10|^(t/))| 2 .) 

IV. Correction of anomalous phase oscillation 

The phase information seems inextricably linked with the 
entanglement information. For the simple case of a fully 
entangled state, either -*=[|00)+e^|ll)] or ^= [|01)+e^|10)], 
our phase indicator returns the cosine squared (of (j) or 6, 



JOURNAL OF JATeX CLASS FILES, VOL. 6, NO. 1, JANUARY 2007 



6 




Fig. 9. The entanglement of 4^[|00> + e^|ll>], as a function of cf), as 
calculated by the QNN, both with and without phase correction as described 
in the text, and compared with Ep (equal to 1 for all </>.) Results are exactly 
similar for the -j= [|01> + e ie 1 10)] states. 



respectively) for the projections, so the procedure is straight- 
forward: Using two copies of a given state, measure the phase 
offset on one copy, then apply a rotation operator to the other, 
and measure the entanglement. Results are shown in Figure 
[9] The linkage does happen in both directions, though. If we 
examine the anomalous oscillation in, e.g., Figure [5] it is also 
remarkably close to a cosine squared. A good approximation 
is shown in Figure [TO] Results are exactly similar for the 
[aoi|01) + e^aiollO)] states. For partially entangled states of 
the type [aoo|00) + aoi|01) + e^an|ll)] (and symmetrical 
equivalents, such as [aoo|00) -f aoi 1 10) + e^aio|10)] ) we also 
get a function of the "contaminant" (for the above state, the 
magnitude aoi) times the cosine squared. This is shown in 
Figure [TT] 

V. Discussion and further work 

We have shown that, for the simple case of a single phase 
offset of a pure state, we can train a quantum neural network to 
correct the anomalous oscillation in an entanglement witness. 
If we know that our experimental setup produces, say, pairs 
of states of the form . 1 [a o|00) + e^a n |ll)] with a 

V a oo+ a n 

high degree of probability, our method would allow for good 
reproducible determination of the entanglement. Of course this 
is not a complete solution even to the measurement of the 
entanglement of a two qubit system. That is, our results in the 
previous section do not transform Figure [7] completely into 
Figure [6| as we have allowed only one of the coefficients at a 
time to be complex. 

We consider this work to be only proof-of-concept: Though 
our phase indicator does give us all relative phases, we 



QNN Entanglement for phase offset Bell states and Curvefitted Model (Approximated) 




Fig. 10. The entanglement of [aoo|00) + e*^an|ll)], as a function of 
4> and aoo> as calculated by the QNN, and compared with a curvefit to 
sin 2 (2aoo) cos 2 (</>) (the partially transparent grid, superimposed.) 



QNN Entanglement of a 00 |00>+a 01 |01>+a 11 e 1 phl |1 1> states 




Fig. 11. The entanglement of [aoo |00)+aoi |01)+e^an 1 11 >] , as a function 
of <f> and aoi, as calculated by the QNN, and compared with a curvefit to 
0.9 cos 2 (1.3aoi) cos 2 (0) (the partially transparent grid, superimposed.) 
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need something more to solve this problem completely, for 
three reasons. First, because each measurement we make on 
a quantum system destroys all other information, so that 
we would need, minimum, four copies of our input state. 
With so many measurements necessary, it seems it might 
be easier simply to determine the state completely and then 
to calculate Ep or whatever we wish to know, analytically. 
Second, inversion of the phase indicator target functions is 
increasingly problematic for the general input state. And third, 
it is not true that if we successfully unwind all the phases 
we will be able to recover the entanglement of the original 
state. Consider the entanglement of formation of the "flat 
state" (equal magnitudes of all of the charge basis states 
| [|00> + e^|01) + e i<9 |10) + e^|ll)]). The concurrence for 
this state is simply C = | sin[ ^~^ +6 ^ ]|, and since Ep is 
a monotonic bijective function of the concurrence, we can 
easily see that (since C is maximum) Ep = 1 whenever the 
difference in phase between the "Bell part" and the "EPR 
part" is 7T. That is, if the "Bell-type" entanglement and the 
"EPR-type" entanglement are exactly in phase, the superpo- 
sition is a product state, and, thus, has zero entanglement; 
if they are exactly out of phase, the density matrix cannot 
be factorized, and, thus, is entangled. This type of phase 
dependence actually causes entanglement, so, if we wish to 
determine the entanglement for a general input state, we will 
have to do more than simply determine all the phases and 
unwind them. It is not yet completely clear how much of 
this behavior the trained QNN net has already captured, since 
for all real coefficients, including negative ones, we get good 
results (Figure [6]). Further work is needed, and is ongoing. 

We have not considered here the multiple qubit systems that 
are of most interest to practical quantum computing; however, 
it is worth considering that a neural approach has natural 
advantages for scale-up. In fact, our entanglement indicator 
did easily generalize from two to three to four to five qubits 
l8lL l9l . We have also not considered mixed states, though 
our entanglement indicator [6] did work also for these types 
of states, which gives us hope that this, larger, work may be 
easily extended to those also. 
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